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Quantum Feedback Control of Atomic Motion in an Optical Cavity 
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We study quantum feedback cooling of atomic motion in an optical cavity. We design a feedback 
algorithm that can cool the atom to the ground state of the optical potential with high efficiency 
despite the nonlinear nature of this problem. An important ingredient is a simplified state-estimation 
algorithm, necessary for a real-time implementation of the feedback loop. We also describe the 
critical role of parity dynamics in the cooling process and present a simple theory that predicts the 
achievable steady-state atomic energies. 
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The control of quantum systems is a problem that lies 
at the heart of several fields, including atom optics and 
nanomechanics. Many of the techniques developed thus 
far for controlling quantum systems in an initially un- 
known state involve the use of dissipative processes to 
achieve a well-defined final state, as in the laser cool- 
ing methods in atom optics pj. A different paradigm 
for cooling involves modifying the system Hamiltonian 
based on information provided by measurements (i.e., 
quantum feedback control) . Such a strategy for control is 
generally applicable to quantum as well as classical sys- 
tems 

0S11HI10I1 As quantum-limited measure- 
ment techniques become more available, quantum feed- 
back control is likely to soon be an essential element in 
atom-optical and nanomechanical toolkits. 

One of the most promising avenues for the study of 
quantum feedback control centers on experiments in cav- 
ity quantum electrodynamics (CQED), where a single 
quantum system may be monitored continuously in real 
time, as has been demonstrated in a series of pioneer- 
ing experiments llfl ITH fl2| . Furthermore, the optical 
potential due to the cavity light provides an easily ad- 
justable means for applying forces on the atom. Feedback 
has been shown to increase the storage time for atoms 
in such an experimental system [l3j |. Despite these ef- 
forts, however, CQED feedback cooling has not yet been 
conclusively demonstrated, highlighting the need for a 
deeper theoretical understanding of this system as well 
as more sophisticated, higher performance cooling algo- 
rithms. Our goal here is to develop approximate estima- 
tion and control algorithms for cooling the atomic mo- 
tion in such a CQED system — an "active cooling" ap- 
proach that is distinct from passive cavity cooling meth- 
ods We will then analyze the performance of these 
algorithms both analytically and through numerical sim- 
ulations. This problem is particularly interesting due to 
the nonlinear form of the measurement operator, a situa- 
tion distinct from previous studies of nonlinear quantum 
feedback control 0, 0|. An important issue that we 
address here is whether the (Gaussian) estimation meth- 



ods developed for linear systems are still effective in this 
problem. 

The setup we are considering is an atom in a micro- 
cavity in the strong-coupling regime, as in the aforemen- 
tioned CQED experiments llfj, llll LLSj , but where the 
output light is monitored via homodyne detection |lfij . 
which gives information about the atomic position in the 
cavity [1?J. We consider the case where the light is reso- 
nant with the cavity, but the detuning from atomic reso- 
nance is much larger than the natural atomic linewidth. 
We thus analyze the coupled atom-cavity system after 
adiabatic elimination of the cavity and internal atomic 
dynamics [18|. The resulting stochastic master equation 
(SME) in Ito form [l9] for the atomic motion conditioned 
on the homodyne measurement, describing the observer's 
best estimate for the atomic state, is 



dp = -i[H eS , p] dt - r[cos 2 (fcX), [cos 2 (fcX), p}] dt 
-y/^{[aa£(kX),p]+ - 2(cos 2 (kX)) p} dW. 



The effective atomic Hamiltonian is 



H e s = 7rP 2 - V max cos 2 (kX), 



(1) 



(2) 



and dW is given in terms of the scaled photodetector 
current (measurement record) df by 



dW = df(t)/^/rj+ y/lkjr (cos 2 (kX))dt. 



(3) 



The physical interpretation of Eqs. and Q is that 
the photocurrent is a direct measurement of cos 2 (kX) 
within the adiabatic approximation. In these expres- 
sions, we are using scaled units such that position is mea- 
sured in units of \Jfif (mw H o)] momentum is measured 
in units of y/Timuj mo , and time is measured in units of 
2ir/ui HO , where u> HO := agky/2h/\mA\ is the oscillation 
frequency in the harmonic approximation for one of the 
optical potential wells, g is the CQED coupling constant, 
m is the atomic mass, A is the detuning from the atomic 
resonance, r\ is the photodetection efficiency, a is the 
mean cavity photon number in the absence of the atom, 
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FIG. 1: Tracking behavior of the Gaussian estimator for the 
position centroid and variance for a sample trajectory. The 
expectation values for the true wave packet (heavy lines) are 
shown with the corresponding Gaussian estimator quantities 
(narrow lines). The cooling algorithm was activated at t = 2. 
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FIG. 2: Evolution of the mean effective energy {H c g) (relative 
to the minimum potential energy), illustrating the cooling- 
algorithm performance. Shown are: cooling via the Gaussian 
estimator (heavy solid), cooling based on the photocurrent 
signal (narrow solid) , and free atomic evolution with no cool- 
ing (dashed). Each curve is an average of 128 trajectories. 



k is the optical spatial frequency, k :— k^/hf (mw HO ), 
Kiax := 7r/fc 2 , T :— 2a 2 g 4 / (A 2 klo ho ) is the effective 
measurement strength, n is the energy decay rate of the 
cavity, and dr(t) = y/ k/oj H o df(t) is the physical pho- 
tocurrent. The model thus has only three independent 
parameters, fc, T, and r/, in addition to any parameters 
that specify how the optical potential is modulated. For 
the atomic parameters, we consider cesium driven on the 
D 2 line, so that m = 2.21xl(T 25 kg and k = 11732 cm" 1 . 
For the cavity parameters, we choose values similar to 
those used in recent CQED experiments pi. Hcl Hl| . yield- 
ing k/2tt = 40 MHz, 5 /2tt = 120 MHz, and a = 1. We 
also assume an optical detuning of A/2tt = 4 GHz to the 
red of the atomic resonance. The corresponding scaled 
parameters are T — 23.6 and k — 0.155. The scaled depth 
of the optical potential is V max = 131, so that the lowest 
two band energies are within 1% of their values in the 
harmonic approximation. We also choose the idealized 
detection efficiency of 77 = 1. However, cooling perfor- 
mance similar to that shown below is achievable even in 
the presence of detector inefficiency (r) well below 0.5) 
as well as other real-world limitations such as hardware 
propagation and computation delays 20] . 

While integration of the SME in order to obtain a state 
estimate is sensible in principle, in practice it is much too 
difficult to integrate the SME in real time. To solve this 
problem we must simplify the estimator. We make a 
Gaussian ansatz for the atomic Wigner function (which, 
equivalently, prescribes how the Weyl-ordered moments 
of the atomic state factorize). In doing so, we reduce 
the effort of integrating a partial differential equation to 
the much more manageable effort of integrating a set of 
five coupled ordinary differential equations for the two 



means (X) and (P) and the three variances V x , V p , and 
C xp := (XP + PX)-(X)(P) H3. The tracking behavior 
of the Gaussian estimator is illustrated in Fig. ^ which 
compares the evolution of the position mean and variance 
for an atom evolving according to Eq. Q with the same 
quantities from a Gaussian estimator evolution beginning 
with a different initial condition. The estimator rapidly 
locks on to the true solution. The ability of the estimator 
to track both quantities well is crucial to the success of 
the control algorithm that we present here. 

Intuitively, an algorithm that responds to the wave- 
packet centroid, raising the optical potential while the 
atom is climbing one side of a well, should cool the 
atom. However, such a strategy eventually breaks down 
due to uncontrolled pumping of energy into the wave- 
packet variances. The key to obtaining an effective 
control algorithm is to consider the energy evolution 
due to the Hamiltonian part of the SME: dt(H c s)ih — 
— {dtV maK )(cos 2 (kX)). From this expression, it is clear 
that a "bang-bang" strategy is optimal for cooling if we 
consider only cyclic modulations of the potential am- 
plitude within a limited range. The potential should 
be switched to the extreme low value when the quan- 
tity — (cos 2 (kX)} is maximized and to the extreme high 
value when this quantity is minimized. We denote these 
extreme values of the potential by (1 — e) 2 V mllx and 
(1 + £) 2 Vmax- Such a control strategy is sensitive to en- 
ergy in both the centroids and variances of the atom. Of 
course, it is not possible to tell whether — (cos 2 (kX)) is 
at an extremum based only on information at the current 
time. Thus, to complete the feedback algorithm, we fit 
a quadratic curve to the history of this quantity at each 
time step and trigger on the slope of this curve at the 
current time. 
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FIG. 3: Evolution of the band populations for the lowest four 
energy bands of the optical lattice, averaged over 128 trajec- 
tories. Top: cooling based on the Gaussian estimator. Bot- 
tom: cooling based on perfect knowledge of the actual wave 
function. 



An example of the cooling algorithm in action is shown 
in Fig. |21 which compares the energy evolution for an en- 
semble of atoms (where each atom is cooled separately 
in an independent simulation) undergoing cooling for 
e = 0.1 with the energy evolution without any cooling 
(i.e., free evolution in the optical potential). Without 
cooling, the atoms simply heat up due to the measure- 
ment back-action (or equivalently, the stochastic cavity 
decay process), until they arc eventually lost from the op- 
tical potential. The heating rate is set by the second term 
in the SME: d t (H cS ) mcas = 87rr(cos 2 (fcA) sin 2 (fcA)). 
The cooling algorithm is able to counteract the measure- 
ment heating, and the atoms settle to a steady-state en- 
ergy much lower than the initial energy. Also shown in 
Fig. |2 is the cooling performance of a much simpler algo- 
rithm, which arises by noting that the photocurrent is a 
direct measure of — (cos 2 (fcA)) (plus noise), and that it 
should be possible to cool the atoms by triggering on the 
photocurrent signal itself (an approach similar in spirit 
to the "differentiating feedback" strategy in Ref. 0). In 
this case, the same quadratic curve-fitting procedure is 
used to temper the measurement noise. Although cool- 
ing to a steady state still occurs in this case, the cooling 
performance is substantially worse than for the Gaus- 
sian estimator. This indicates that the Gaussian esti- 



mator acts as a nearly optimal filter (in the same sense 
as a Kalman filter \22\) for the relevant information in 
the photodetector signal. The results shown were gener- 
ated assuming initial conditions and curve-fit details as 
follows, although the results do not sensitively depend 
on these values: The atoms begin in a coherent state 
{V x = V p = 1/2), and each has an initial centroid energy 
corresponding to a displacement from the well center of 
58% of the distance to the edge of the well, chosen so 
the atoms are clearly trapped but not particularly cold. 
The initial locations of the atoms are distributed uni- 
formly between the bottom of the well and this maxi- 
mum distance. The Gaussian estimator begins with the 
(impure-state) initial conditions (X) c — 6, (P) c = 0, 
V x = Vp = 1/V2, and C c — 0. The estimator is evolved 
in steps of At = 0.0005, and at each step the quadratic 
curve is fitted to the values of — (cos 2 (fcA)) computed 
from the 300 most recent estimator steps. 

The cooling performance is further illustrated in Fig.|3 
which shows the populations in the lowest four bands of 
the optical potential. Most of the steady-state popula- 
tion is in the lowest two bands, accounting for 94% of 
the population, and these two bands are almost equally 
populated. We now show that this behavior can be un- 
derstood in terms of the parity of the atomic state. The 
atomic parity is invariant under Hamiltonian evolution 
in the atomic potential (and hence the influence of the 
control algorithm). However, the measurement term in 
the SME causes the atomic parity to diffuse according to 

d(V) = -y^8rjr \(Vcos(2kX)) - (V)(cos(2kX))] dW, 

(4) 

where V is the parity operator. In fact, this is an unbi- 
ased, nonstationary diffusion process that causes (P) to 
"purify" at late times to the extreme values ±1 of pure 
parity, so that the homodyne measurement acts like a 
QND measurement of the atomic parity. Since the atomic 
initial condition described above corresponds to (V) = 0, 
the measurement process drives the atoms to states with 
either pure even or pure odd parity with equal probabil- 
ity. The obvious consequence of this effect for cooling is 
that only half of the atoms can be cooled to the ground 
band, while the remaining half can be cooled at best to 
the first excited band [2l|. However, due to this purifica- 
tion effect (which means that the observer knows the final 
state of the atom on each individual run) , cooling to the 
first excited band is essentially equivalent to cooling to 
the ground state, since the atom can be driven between 
these two bands with an additional coherent process (e.g., 
modulation via an additional classical optical potential). 

Much of the residual energy (excitation out of the low- 
est two bands) is due to the inability of the Gaussian 
estimator to perfectly track the true atomic state. In 
this light, the cooling performance is quite impressive, 
especially considering the highly non-Gaussian nature of 
the first excited state, toward which half the atoms are 
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FIG. 4: Final energies, measured over the interval from t = 90 
to 100, as the switching amplitude e varies. Circles: cooling 
based on the Gaussian estimator. Squares: cooling based on 
perfect knowledge of the actual wave function. Dashed line: 
simple cooling theory, Eq. (J^J. Inset: magnified view of the 
same data. Error bars reflect standard errors from averages 
over 128 trajectories. 



cooling. To show this behavior explicitly, we can remove 
the estimator from the problem, and repeat the same 
ensemble of simulations where the cooling algorithm op- 
erates based on the value of — (cos 2 (kX)} computed with 
respect to the actual atomic wave function. These results 
are also shown in Fig. |21 and in this case the lowest two 
bands account for 98% of the atomic population. 

The steady-state atomic energies can be well described 
by a simple analytic theory, which we now derive. We will 
assume that the atomic energy is near the ground-state 
energy Eq, so that we can make the harmonic approxi- 
mation for the potential. We will further assume for the 
moment that all the energy AE in excess of Eq is associ- 
ated with the wave-packet centroid (i.e., the wave packet 
is not squeezed). Then the quantity — V max (cos 2 (A:X)) 
exhibits oscillations of peak-to-peak magnitude AE. As- 
suming that the switching occurs at the extremal times, 
the energy-evolution expression above leads to a coarse- 
grained cooling rate of <9 t (-ff e ff)fb = — 8e((i? ff) — ^o)- 
On the other hand, the heating rate discussed above can 
be rewritten as <9 t (i/ c ff) mcas = 4Tk 4 (H c g). Steady state 
occurs when these two rates sum to zero: 



Eq + Ei 

2(1-/3)' 



(5) 



Here, (3 :— Tk 4 /2e, and we have made the replacement 
Eq — > (Eq + Ei)/ 2 in view of the above parity consider- 
ations, where E\ is the energy of the first excited band. 

We can also derive a theory for the other ex- 
treme, where we assume all the excess energy is as- 
sociated purely with squeezing of the wave packet 
and not with the centroid. In this case, the oscilla- 
ti ons of — V max ( cos 2 (kX)} are of peak-to-peak magnitude 
^(iJcff) 2 - E$. Then the (l-/?)" 1 dependence in Eq. © 
is replaced by (1 



this theory should be valid, this prediction is always be- 
low that of Eq. J5J, and thus the centroid-only theory 
more faithfully represents the physical cooling limit. 

Eq. J5| is compared with simulated late-time energies 
in Fig. 01 The simulated energies are in good agree- 
ment with the simple theory, although they are gener- 
ally higher than the simple theoretical prediction. Again, 
much of this discrepancy is due to the imperfect tracking 
behavior of the Gaussian estimator, and the simulated 
energies from cooling based on perfect knowledge of the 
actual wave function are in better agreement. Qualita- 
tively, we see a transition between ineffective and effective 
cooling with e; essentially, e must be large enough that 
the cooling rate exceeds the heating rate, ensuring that 
the system is controllable. Beyond this transition, the 
system displays a remarkable insensitivity to the value of 
the switching amplitude; we observe similar robustness 
to the changes in the other parameters T, k, and 77. 
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